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In order to single out dominant phenomena that account for carrier-controlled magnetism in p- 
Cdi_2,Mn2,Te quantum wells we have carried out magneto-optical measurements and Monte Carlo 
simulations of time dependent magnetization. The experimental results show that magnetization 
relaxation is faster than 20 ns in the paramagnetic state. Decreasing temperature below the Curie 
temperature Tc results in an increase of the relaxation time but to less than 10 /iS. This fast 
relaxation may explain why the spontaneous spin splitting of electronic states is not accompanied 
by the presence of non-zero macroscopic magnetization below Tc- Our Monte Carlo results repro- 
duce the relative change of the relaxation time on decreasing temperature. At the same time, the 
numerical calculations demonstrate that antiferromagnetic spin-spin interactions, which compete 
with the hole-mediated long-range ferromagnetic coupling, play an important role in magnetization 
relaxation of the system. We find, in particular, that magnetization dynamics is largely accelerated 
by the presence of antiferromagnetic couplings to the Mn spins located outside the region, where 
the holes reside. This suggests that macroscopic spontaneous magnetization should be observable 
if the thickness of the layer containing localized spins will be smaller than the extension of the hole 
wave function. Furthermore, we study how a spin-independent part of the Mn potential affects Tc. 
Our findings show that the alloy disorder potential tends to reduce Tc, the effect being particularly 
strong for the attractive potential that leads to hole localization. 

PACS numbers: 75.50.Pp,05.10.Ln,75.30.Et,78.55.Et 



I. INTRODUCTION 

Modulation-doped p-type (Cd,Mn)Te/(Cd,Mg,Zn)Te 
quantum wells (QWs) remain a unique medium allow- 
ing to probe carrier-induced Ising-like ferromagnetism in 
the two-dimensional (2D) case,^ 2345 gg j^-^ ^j^jg system 
the mean free path is longer than the QW width. The 
presence of ferromagnetism was revealed by the obser- 
vation of spontaneous splitting of the photoluminescence 
(PL) spectrum, which indicated that local spin order- 
ing is stable, below the Curie te mp erature Tc, for times 
longer than the exciton lifetime.'^El Theoretical analyzes 
have shown that the magnitude of Tq as a function of 
the Mn concentration x and the hole areal density p can 
be quantitatively described provided that in addition to 
carrier-mediated ferromagnetic spin-spin couplings, the 
presence of competing short-range antiferromagnetic su- 
perexchange interactions as well as hole correlation ef- 
fects are taken into account .'^^^'^'^ These experiments 
and their quantitative interpret ations were carried out 
prior to later theoretical studies. ^^^1111111 

A surprising result of both optical and magnetic mea- 
surements is the absence of hysteresis loops and, he nce, 
of macroscopic spontaneous magnetization below Tc-'^ 
This finding questions our understanding of the actual 
ground state of the system. It has been suggested, in par- 
ticular, that the formation of spin density waves, driven 
by the g-dependent carrier susceptibility, may account for 



the absence of spontaneous magnetization .^^l This possi- 
bility has recently been considered in the context of the- 
oretical search for carrier-mediated nuclear magnetism in 
2D systems .ISl 

In order to determine the importance of various 
phenomena that control magnetism in such reduced- 
dimensionality magnetically disordered systems we have 
carried out time-resolved magneto-optical measurements. 
Furthermore, with the goal to obtain information on 
mechanisms controlling spin dynamics, we have extended 
our previouPI Monte Carlo (MC) simulations of magneti- 
zation in p — Cdi_a;Mna;Te QW. The experimental find- 
ings provide a dependence of magnetization relaxation 
on temperature and carrier density. Above Tq the relax- 
ation time is found to be shorter than 20 ns, while below 
Tc magnetization persists up to a few microseconds, a 
time scale consistent with the absence of spontaneous 
magnetization in the static measurements. Interestingly, 
we have obtained a similar relative prolongation of the 
relaxation time below Tq from MC simulations based on 
the Metropolis algorithm and on the determination of 
the one-particle hole eigenfunctions at each MC sweep. 
The simulations explain also the rate of increase of the 
relaxation time with the carrier density. We conclude 
that the absence of magnetic hysteresis in this 2D Ising 
system can be explained without taking carrier-carrier 
correlation into account. 

At the same time, by analyzing results of our Monte 
Carlo simulations, we are able to show that short range 
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spin-spin antiferromagnetic (AFM) interactions play a 
crucial role in accelerating magnetization dynamics. This 
corroborates the outcome of the previous results showing 
a narrowing of the hysteresis loop by AFM interactions P 
However, we realize now, based on much more exten- 
sive simulations, that the effect of AFM interactions be- 
comes much reduced if the thickness of the layer contain- 
ing Mn spins is taken to be smaller than the extent of 
the hole wave function. This finding, from the one hand, 
substantiates theoretical considerations of Boudinet and 
Bastard,!^ suggesting that magnetization relaxation of 
bound magnetic polarons in p-type CdTe/(Cd,Mn)Te 
QWs occurs owing to the AFM coupling to the Mn spins 
located outside the relevant Bohr radius. On the other, 
this result implies that macroscopic spontaneous magne- 
tization and magnetic hysteresis should be recovered in 
quantum wells, in which the thickness of the Mn layer 
would be smaller than the region visited by the holes. 

Furthermore, we study how a spin-independent part 
of the potential introduced by Mn impurities affects Tq ■ 
Our simulations show that alloy disorder tends to reduce 
Tc. The effect is particularly dramatic for the attrac- 
tive alloy potential which, if sufficiently large, leads to 
a strong hole localization. This result substantiates the 
notion that delocalized or weakly localized carriers are 
necessary to generate a sizable ferromagnetic coupling 
between diluted localized spins. 

This article is organized as follows. In Sec. [ll] we 
present experimental findings of magneto-optical studies 
on magnetization dynamics, which have motivated the 
theoretical effort. The theoretical model that describes 
the carrier-induced ferromagnetism in diluted magnetic 
semiconductor (DMS) quantum wells is introduced in 
Sec. In Sec. |IV| we present MC investigations of the 
dependence of the critical temperature on the carrier den- 
sity and spin-independent alloy potential. Section [V] con- 
tains results of our MC studies of spin dynamics and 
of spin ordering, which serve to interpret the findings 
of magneto-optical measurements. Specifically, subsec- 



tion VA is devoted to magnetization dynamics, while 
subsection |VB| contains a discussion of possible domain 
sizes. Our conclusions are summarized in Sec. IVII 



II. EXPERIMENTAL RESULTS 

The experimental results have been obtained for sam- 
ples grown in Grenoble by molecular beam epitaxy on 
Cdn.88Zno.12Te (001) substrates.^ The modulation- 
doped p-type structures contain a single 8 nm wide 
Cdi_:rMn:rTe QW embedded by Cdi^y^MgyT^n^Te bar- 
riers, in which the Mg content {y = 0.25-0.28) results in a 
sizable valence band offset, while the presence of Zn (z — 
0.07-0.08) ensures a good lattice match to the substrate. 
The front barrier is doped by nitrogen acceptors. The 
distance between the QW and the doping layer is 20 nm, 
which results in the hole density up to 3 x lO^-'^cm"^.'^^ 
These holes occupy the ground state subband whose scat- 



tering broadening is much smaller than the distance to 
the next subband, making the system to be truly 2D. 

According to the well estabhshed procedure, 
magnetic properties of the QW are probed by 
magnetospectroscopy.l^J^ The photoluminescence 
(PL) is excited by a cw laser with the photon energies 
below the absorption edge of the barrier material. A typ- 
ical PL spectrum of a paramagnetic p-type QW consists 
of a single line related to the radiative recombination of 
charged excitons (X+). Depending on its polarization 
chirality, the emission line exhibits a downward or 
upward shift in the magnetic field. Owing to a strong 
sp-d exchange interaction specific to DMSs, this Zeeman 
shift is giant (typically over 20 meV in 1 T), and its 
magnitude is proportional to the local magnetization 
of the Mn spins. Importantly, a combination of strain, 
confinement, and spin-orbit interaction makes that the 
giant shift vanishes virtually entirely for the in-plane 
magnetic field. Below a critical temperature, which we 
identify as the Curie temperature Tq, a spontaneous 
splitting of the PL line into two components is observed 
in the absence of an external magnetic field. The energy 
distance between these two lines, or the downward shift 
of the lower line, provides information on an average 
value of spontaneous magnetization in the regions visited 
by the holes at given temperature. 

Surprisingly, the emitted light at the wavelength corre- 
sponding to either of these two lines remains unpolarized, 
even after ramping the field down to zero below Tq- A 
magnetic field of the order of 20 mT, much higher than 
the demagnetization field for this diluted alloy, is nec- 
essary to achieve a full circular polarization of the lines 
at T « O.STc.l^This indicates the absence of a macro- 
scopic spontaneous magnetization at these low tempera- 
tures, the conclusion confirmed by a direct magnetization 
measurements.'^ Moreover, even focusing detection on 
the spot with diameter below 1 /im does not result in the 
circularly polarized emission. This could be explained ei- 
ther by a small size of the relevant magnetic domains or 
by fast fluctuations of the magnetization direction. Al- 
ternatively, a spin-density wave could be considered as 
the relevant ground state of this reduced dimensionality 
correlated system. However, a selective excitation of PL 
by circularly polarized light leads to circularly polarized 
PL.I^l This indicates that local spin ordering is stable for 
times longer than the exciton lifetime, estimated to be in 
the 100 ps range. 

Some of us have recently developed a technique for 
studies of magnetization dynamics by means of time re- 
solved PL measurements performed after a short pulse 
of the magnetic field.'* ^ Pulses of the magnetic field are 
produced by a magnetic coil mounted at the surface of 
the sample. The illumination with a laser beam and the 
collection of the PL signal are performed along the axis 
of the coil (Faraday configuration) . The coil diameter of 
0.5 mm results in a small value of inductance and allows 
to obtain short rise and fall times of about 10 ns. A 
2 A current produces a magnetic field of about 40 mT. 
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FIG. 1: Experimental: Temporal decay of magnetization in 
the Cdo.96Mno.04Te quantum well at various temperatures 
determined as a difference in photoluminescence intensities 
collected for two circular polarizations after a pulse of the 
magnetic field. The Curie temperature is 2.5 K. 
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FIG. 2: Experimental: Temporal decay of magnetization in 
the Cdo.96Mno.04Te quantum well determined as a difference 
in photoluminescence intensities for two circular polarizations 
collected at various wavelengths after a pulse of the magnetic 
field. The Curie temperature is 2.5 K. 



The time evolution of the PL during and after the pulse 
is probed with resolution down to 10 ns. Magnetization 
dynamics is studied by analyzing the difference in PL in- 
tensities for two circular polarizations after the pulse of 
the magnetic field. Results of measurements shown here 
have been obtained at different temperatures below and 
above Tc, which is 2.5 K for the studied Cdo.96Mno.04Te 
QW. 

Above Tc, in the paramagnetic state, a single PL line 
is observed at zero field, and the signal induced by the 
magnetic pulse is related to Zeeman splitting of this line. 
Its relaxation time is found to be shorter than 20 ns. Such 
a value is just t o be expected for (Cd,Mn)Te containing 
4% of Mn.™ 

Typical temporal profiles for three different tempera- 
tures below Tc — 2.5 K are presented in Fig. [l] The 
relaxation time becomes much slower when decreasing 
temperature (see Fig. [l]). Magnetization persists up to a 
few microseconds when T « 0.5Tc. However even for the 
lowest temperatures, the relaxation time is shorter than 
10 IIS, and under all experimental conditions the rema- 
nent magnetization vanishes totally after 20 fis. This 
time scale is much shorter than an acquisition time of 
standard optical measurements, which accounts for un- 
measurably small values of the coercive field observed 
experimentally. 

It is interesting to note that the decays are not mono- 
exponential. The relaxation begins with a rapid decay, 
faster than 100 ns, but it also contains a slower compo- 
nent, in the microsecond range. One could expect that 
regions of higher hole density have higher local magne- 
tization resulting in higher barriers for the magnetiza- 
tion reversal. In order to examine this possibility we 
carried out relaxation measurements for a constant tem- 



perature (below Tc), and different detection wavelength 
which was tuned over the low energy component of the 
PL line. As mentioned above, this line shows a red shift 
when decreasing temperature below Tc- This shift is 
proportional to Mn magnetization, who se m agnitude is 
primarily determined by the hole density.'^'i^ Due to spa- 
tial fluctuations in the hole density,'^ the line is signifi- 
cantly broadened. The low energy side of the line comes 
from the regions containing higher hole densities, while 
the high energy side corresponds to locations visited by 
fewer holes. This provides an opportunity to trace the 
relaxation time as a function of the hole concentration. 
The obtained temporal magnetization profiles are pre- 
sented in Fig. [2] It is seen that the contribution of the 
slow component increases for the regions containing a 
high hole density. This finding as well as the observed 
temperature variation of the magnetization decay time 
will be compared to results of numerical simulations in 
the next section. 



III. THEORY AND SIMULATION PROCEDURE 

To understand the experimental results described in 
Sec. |TT] we carried out Monte Carlo simulations for the 
subsystem of Mn spins which are at the same time subject 
to short-range antiferromagnetic and long-range hole- 
mediated ferromagnetic spin-spin interactions. Accord- 
ing to our results, the presence of these two kinds of com- 
peting interactions accounts for a non-standard character 
of ferromagnetism in the studied system. Since we are 
mainly interested in modulation-doped p-Cdi_a;Mna;Te 
quantum wells, all our results are obtained using the 
parameter values characteristic for this system, but the 
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methodology we use may as well be applied to other 2D 
p-type II- VI DMS structures. 

The antiferromagnetic interactions originate from su- 
perexchange, the mechanism specific to undoped charge 
transfer insulators, such as (Cd,Mn)Te. The character 
and magnitude of these interactions are known from pre- 
vious extensive experimental studies of magnetic proper- 
ties of (Cd,Mn)Te. In order to find out how the presence 
of holes affects the system we compute, in the spirit of the 
adiabatic approximation, the total energy of the holes at 
a given configuration of the Mn spins. In this way our 
approach encompasses automatically the description of 
carrier-mediated exchange interactions within either p-d 
Zener model or Ruderman-Kittel-Kasuya-Yosida theory. 
The carrier total energy is evaluated neglecting hole cor- 
relations, the approximation that leads to an underesti- 
mation of the ferromagnetic coupling.l^ Accordingly, our 
computed Curie temperatures are expected to be system- 
atically lower than those observed experimentally. Fur- 
thermore, in accord with the effective mass theory, and 
owing to a relatively small disorder in our structures, 
we assume that the in-plane and the perpendicular hole 
motions can be factorized, i.e., that the total envelope 
function assumes the form \E'o-(i?) = 'ijj^(r)ip{z), where 
R = (r, z) with in-plane and perpendicular positions de- 
scribed by r = (x, y) and z respectively; ^{z) is the wave 
function of a hole residing in the ground state subband of 
the QW in question. For the biaxial strain under consid- 
eration, the relevant subband is of the heavy hole origin, 
so that the spin-orbit interaction fixes the hole spin in 
the direction of the growth axis. Furthermore, since the 
heavy holes penetrate only weakly the barriers, '^{z) in 
the form corresponding to an infinite QW is adequate for 
our computations. 

Within the above model and in the absence of an ex- 
ternal magnetic field, the Hamiltonian Ti. of our system 
takes the form, 

^ = ^ + '^pd + '^dd- (1) 

This Hamiltonian describes noninteracting carriers in the 
heavy hole parabolic valence subband for which the in- 
plane effective mass is m* — 0.25mo for the QW in 
question .^12] Since the coupling between the holes and the 
Mn ions is weak in (Cd,Mn)Te,'22l we can consider the p-d 
interaction in a simple contact form, 

n^d = \P]^S,_,5{r - n) , (2) 

i 

where the hole spin j = 3/2, = ±3/2, is coupled to the 
Mn ions randomly distributed over the fee lattice sites r^. 
We treat Mn spins Si as classical vectors with S — 5/2. 
While this approximation is certainly qualitatively valid 
for the large spin in question, quantitatively - according 
to the mean-field theory - it leads to an underestima- 
tion of the Tq value by {S + ^)/ S. For the p-d exchange 



integral we take j3 = —5.96 x 10~^ eV nm^, which corre- 
sponds to the exchange energy PNq = —0.88 eV, where 
A'o is the cation concentration.l^il 

We have also tested how our results are affected by 
the spin-independent part of the potentials introduced 
by Mn ions. To take this alloy disorder into account we 
add to Eq. |2] a term 

nAB = 5Vj2Kr-n)\ip{z,)\^ (3) 

i 

with 6V ~ 3.93 x lO^'^ eV nm'^ which corresponds to 
the valence band offset WNq = -SVNq = -0.58 eV in 
(Cd,Mn)Te.22 23 jt is known that the effect of the Mn 
isoelectronic impurities on the carrier wave function is 
governed by the ratio of the total Mn potential U to its 
critical value C/c < at which a bound state starts to 
form,l2ni 

U/Uc = 6m*[W - {S+l)l3/2]/{7r^h^b), (4) 

where b is the potential radius. In the case under consid- 
eration, the hole spin is fixed along z-axis, and Mn po- 
larization (Sz) is typically below 10%. Accordingly, for 
S'-l- 1 — > {Sz) the total Mn potential is expected to be ef- 
fectively repulsive for the holes in (Cd,Mn)Te, U/Uc < 0. 
Interestingly, in many important systems, such as oxides 
and nitrides, the spin-independent part of the magnetic 
impurity potential is actually attractive.^^Sl Accordingly, 
we have also performed some simulations of magnetiza- 
tion for SV = —3.93 x 10~^ eV nm'^ which corresponds 
to the valence band offset WNq = -SVNq = -^0.58 eV. 
However, it should be noted that by using the (5-like po- 
tentials, which corresponds to 6 0, we actually overes- 
timate U/Uc and determine an upper limit of the effect 
of the alloy potential upon Tq. With this in mind, we 
have carried out all MC dynamics studies presented in 
Sec. El for SV = 0. 

The short-range intrinsic antiferromagnetic (AFM) in- 
teraction between the Mn spins is given by, 

Wdrf = —2^5 ^ JijSi ■ Sj, (5) 

ij 

with Jij = —6.3, —1.9, —0.4 K for the nearest, next near- 
est, and next next nearest neighbors, respectively.^^S Our 
MC simulations confirm that these values of Jij repro- 
duce correctly the temperature dependence of the spin 
susceptibility in the absence of holes but we note that 
other sets of Jij values have also been proposed in the 
literature.l^ 

The above approach disregards the presence of long- 
range dipole-dipole interactions between Mn spins. Both 
theoretical evaluations and MC simulations we have car- 
ried out show that the dipole-diploe interactions can be 
safely ignored in the case under consideration. 

In order to determine hole eigenfunctions and eigen- 
values for a given configuration of Mn spins we assume 
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periodic boundary conditions in the QW atomic lay- 
ers and diagonalize 7i in a plane-wave basis with 2D 
wave vectors truncated at a radius — 5.52(271/ L) or 
kc = 8.10(27r/L)P This part of the 2D fc-space cor- 
responding to 97 or 213 fc-states, respectively, is suffi- 
cient to ensure the convergence for all hole density values 
considered here. The energy of the holes is determined 
by summing up the lowest eigenvalues corresponding to 
a given number of holes N^. This procedure assumes 
that the hole gas is degenerate, i.e., neglects thermal 
broadening of the hole distribution, the effect considered 
previously,^ and found to be small. In order to keep fc- 
space shells filled up, we have considered = 1, 5, 9, 
13, 21, 25, 29, 37, 45, 49, and 57, which correspond to 
values of hole concentrations p up to 1.11 x 10^"'^ cm~^ 
for the chosen size of the simulation box. 



In our MC simulations we employ the Metropolis al- 
gorithm in which the hole eigenfunctions and eigenvalues 
are updated at each MC sweep (in one MC sweep all Mn 
spins are rotated). This procedure if applied after each 
single spin rotation would be computationally excessively 
time consuming and, therefore, we followed the idea of 
the so-called perturbative Monte Carlo method. asHSl 
We typically keep 2000 initial MC sweeps (to thermalize 
the system), followed by lO^'-lO^ MC sweeps used for a 
further analysis. Most of the numerical results presented 
in this paper have been obtained for the magnetic field 
put to zero. However, our software allows as well simula- 
tions in non-zero magnetic fields. We have occasionally 
used this opportunity to obtain some of the initial con- 
figurations. 



For the computations presented here we adopt the sim- 
ulation box of dimensions L x L x Lw , where L < 350ao 
and Lw = Soo, with the fee lattice constant Oq = 
0.647 nm in our case. Hence, the quantum well width 
is Lw — 5.2 nm and the area L^ =(226 nm)^. We 
note that a typical QW width of experimental samples 
is somewhat grater, L"^^ = 8 nm. However, the adop- 
tion of a smaller value of Lw allows us to treat systems 
with a larger area and, at the same time, since within 
the mean-field theorjJ^ cx 1/Lw, results in a partial 
compensation of a systematic error stemming from treat- 
ing the spins classically, L'^^S/[Lw{S + 1)] = 1.1. We 
randomly occupy 4% of cations sites, which corresponds 
to 156800 Mn spins. In view of the required accuracy, 
this large system size makes averaging over the Mn dis- 
tribution and extrapolation to an infinite system size un- 
necessary. The number of Monte Carlo sweeps needed 
for the convergence is of the order of 10^ to 10^. Hence, 
we have to use a (pseudo) random number generator with 
a sufficiently long period. After testing various genera- 
tors, the Mersenne Twister developed by Matsumoto and 
Nishimura has been selected, because of its huge period 
and its very high order of dimensional equidistribution.'^'^ 



IV. CURIE TEMPERATURE 

The critical temperature Tc has been determined in 
the standard way: by MC simulations we have calculated 
the temperature dependence of the spin projections and 
the spin susceptibility for both hole and magnetic ion 
spins. To denote spins we use the notation a = s (for 
holes) and a — S (for Mn ions). Similarly, N — (for 
the number of holes) and N = Nun (for the number of 
Mn ions). At each MC sweep, the spin projections per 
hole or per Mn ion are given by: 



N ' 



J = x,y,z. 



(6) 



We symbolize statistical averages of the above quantities 
by (...).'^In other words, e.g. 



nsum 



and 



(Kl> = 



yr. 



(7) 



(8) 



where n denotes successive MC sweeps used for the statis- 
tical average and nsum denotes the total number of MC 
sweeps used for the statistical average. The spin suscep- 
tibilities per spin are defined in the following manner.*^ 



1 



(9) 



We associate the temperature at which the magnitude 
of spin susceptibility reaches a maximum with the Curie 
temperature Tc of the system. Usually, the MC com- 
putations have been performed with the step of 0.2 K. 
Therefore, an uncertainty in Tq values is usually at least 
±0.1 K. Some examples of the temperature dependence 
of the spin projections and of the spin susceptibilities are 
presented in Fig. |3] Specifically, we have plotted {\Sz\) 
and xs for the Mn ions (Figs. [3ja) and ^c)) as well 
as (|sh|) and Xh for the holes (Fig. ^h) and[3|d)), for 
the hole numbers iVh = 13 and 29, which correspond to 
p = 2.54 X 10"^° and 5.66 x 10^" cm~^, respectively. In 
Fig. |3] the case when the alloy disorder is absent (no AD) 
is compared to the situation of the repulsive (-fAD) al- 
loy disorder potential. In Fig. |4] we show the effect of 
the attractive (-AD) alloy disorder potential by referring 
to the data obtained without alloy disorder (noAD), for 
Nh = 13 and 45, i.e., p = 2.54 x 10^° and 8.78 x lO^" 
cm^^, respectively. 

In Fig. |5] we summarize our MC results for Tc as a 
function of the hole concentration obtained neglecting 
alloy disorder (SV = 0) as well as for SV — ±3.93 x 
10~^ eV nm^, corresponding to the repulsive and at- 
tractive alloy potential introduced by Mn, respectively. 
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FIG. 3: Monte Carlo simulations: the z-component of the absolute spin projections {\Sz\) and {|sz|) (open symbols - dotted 
lines; left scale) and the spin susceptibilities xs and Xs (full symbols - solid lines; right scale) for the Mn ions (a,c) and for the 
holes (b,d), in the absence of alloy disorder (noAD, SV — 0), and for the repulsive (+AD, 5V > 0) alloy disorder potential. A 
maximum in the temperature dependence of spin susceptibility is identified as the Curie temperature below which the holes 
induce a long range ferromagnetic order. Here the number of holes A^h = 13 (a,b) and 29 (c,d), i.e., p = 2.54 x 10^'^ and 
5.66 X 10^° cm~^, respectively. 



We see that the attractive potential washes the ferro- 
magnetism virtually entirely out. A contour plot of hole 
eigenfunctions demonstrates that holes get localized for 
the chosen magnitude of 6V. Our results substantiate 
therefore the view that delocalized or weakly localized 
carriers are indispensable to set a long-range order be- 
tween diluted spins. Obviously, smaller amplitude of the 
attractive potential or greater values of hole concentra- 
tions will lead to carrier delocalization and the reentrance 
of a ferromagnetic order. 

Interestingly, as seen in Fig. [5] the repulsive potential 
SV > 0, the case of (Cd,Mn)Te, also leads to reduced 
magnitudes of Tq comparing to the values determined 
for SV — 0. We interpret this finding by noting that in 
the presence of a repulsive potential, the amplitude of 
the wave function at Mn ions is diminished comparing to 
the case SV = 0. This effectively reduces the p-d cou- 
pling and shifts the appearance of the carrier-mediated 



ferromagnetic order to lower temperatures. 

It is instructive to compare these MC findings 
to the expec tations of the mean-field approximation 
(MFA),1EEE1 derived neglecting entirely chemical and 
thermal fluctuations. We see in Fig. |5] that, somewhat 
fortunately, the MC data for SV = and MFA results 
for classical spins (solid line) are in a good quantitative 
agreement at high hole densities. In particular, both ap- 
proaches predict that the magnitude of Tq does not vary 
with the hole density, the result reflecting the energy in- 
dependence of the density of states in the 2D case. How- 
ever, at low hole densities, Tc values obtained from the 
MC simulations tend to decrease, the effect associated 
with a broadening of the density of states induced by 
hole scattering on Mn spins® ^ and encompassed by our 
MC simulations. The scattering-induced lowering of Tq 
is even more apparent if SV > 0. 

Owing to mutual compensations between ferromag- 
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FIG. 4: Monte Carlo simulations: the z-component of the absolute spin projections (l^zl) and {\sz\) (open symbols - dotted 
lines; left scale) and the spin susceptibilities xs and Xs (full symbols - solid lines; right scale) for the Mn ions (a,c) and for the 
holes (b,d), in the absence of alloy disorder (noAD, 5V — 0), and for the attractive (-AD, SV < 0) alloy disorder potential. A 
maximum in the temperature dependence of spin susceptibility is identified as the Curie temperature below which the holes 
induce a long range ferromagnetic order. Here the number of holes A^h = 13 (a,b) and 45 (c,d), i.e., p = 2.54 x 10^" and 
8.78 X 10^° cm"^, respectively. 



netic and antiferromagnetic interactions, the magnitude 
of Tc is rather sensitive to the presence of hole corre- 
lations. This is shown by dashed line in Fig. [5] which 
has been obtained by incorporating into theory within 
the MFA the Fermi liquid Landau parameter Ap — 2. 
This procedure increases the ferromagnetic contribution 
to Tc by a factor of two but since carrier correlations do 
not affect the compensating AFM term, the resulting in- 
crease of Tc is actually much greater. Therefore, in order 
to compare MC results to experimental findings we treat 
Ap as a fitting parameter. The theoretical values of Tc 
displayed in Fig. [g] has been determined from; ^-^-'''^' 

x{Tc)=x{t[^''^)/Af. (10) 

Here x(T) ^ T~" is the spin susceptibility of the an- 
tiferromagnetically coupled Mn spins in the absence of 
the holes, where a = 0.79 according to our MC simula- 



tion, and is the Curie temperature at a given hole 

concentration, determined from the MC simulations that 
neglect carrier correlations {Ap — 1). We see that for 
p — I X 10^^ cm^^ by taking Ap — 1.5 and Ap — 3 for 
the case 5V = and 6V — 3.93 x 10^^ eV nm'^, respec- 
tively we obtain a good description of the experimental 
findings. In particular, a decrease of Tc values with low- 
ering of the hole density observed experimentally^ and 
also visible in our MC simulations, results from scatter- 
ing broadening of the density of stateSjf^ particularly rel- 
evant at low Fermi energies. However, the experiment 
implies that this lowering extents to higher concentra- 
tions than the range predicted by the simulations, which 
indicates that additional scattering mechanisms operate 
in real samples. 

Having discussed the ferromagnetic ordering tempera- 
ture Tc , we turn to the temperature dependence of Mn 
and hole spin projections, as depicted in Figs. [3] and [4] 
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FIG. 5: Curie temperature Tc as a function of the hole area 
density p. Sohd symbols show results of Monte Carlo simula- 
tions carried out for classical spins, QW width Lw ~ 5.2 nm, 
and Mn content 4% neglecting carrier correlations in the ab- 
sence of alloy disorder {SV — 0, squares) as well as for re- 
pulsive and attractive alloy potentials {SV > 0, triangles 
and SV < circles, respectively). Results obtained within 
the mean- field approximation (i.e., disregarding thermal and 
chemical fluctuations) are shown by the solid line. The dashed 
line presents the MFA results taking into account carrier cor- 
relations with the Landau parameter Ap — 2. 



We see that the holes become entirely spin polarized, 
(|sz|) ^ 1/2, immediately below Tq- On the other hand, 
the increase of the Mn spin projection ( 15*2 1) on lower- 
ing temperature is much slower. This is because, the 
molecular field produced by the spin polarized carriers is 
typically below 1 kOe, much too small, even at 0.1 K, 
to saturate Mn spins which are coupled by AFM inter- 
actions. 

In Fig. [T] we compare the magnitudes of the Mn spin 
projection {\Sz\) from our MC simulations at 0.1 K with 
the values expected from the MFA for the saturated car- 



rier spins,' 



ilSMT) 



(11) 



where x(2^) ^"^d Xo(T) are the values of Mn spin sus- 
ceptibilities per Mn ion computed without holes, in the 
presence and in the absence of the AF interactions, re- 
spectively. We see that the MFA quite correctly repro- 
duces the computed magnitudes of {\Sz\) as a function of 
the areal hole concentration p. Furthermore, while the 
values of Tq are systematically smaller in the presence 
of the repulsive alloy potential comparing to the case 
SVNq = 0, the magnitude of {\Sz\) is seen to be larger 
in this case. We assign this surprising result to the fact 
that the repulsive alloy potential delocalizes majority- 
spin holes, which can therefore embrace more effectively 



FIG. 6: Curie temperature Tc as a function of the hole areal 
density p. Solid symbols show results of Monte Carlo simula- 
tions carried out for classical spins, QW width Lw ~ 5.2 nm, 
and the Mn content 4% in the absence of alloy disorder 
{SV = 0, squares) and for repulsive alloy potentials {SV > 0, 
triangles) obtained with the Fermi liquid Landau parameters 
Ap = 1.5 and Ap = 3. Experimental resultJ^ for QW width 
Lw ~ 8 nm and Mn content 3-5% for two samples (open 
circles and triangles). 
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FIG. 7: Monte Carlo simulations: z-component of the abso- 
lute spin projection of Mn spins (ISzl) at 0.1 K in the absence 
of alloy disorder {SV — 0, squares) and for repulsive alloy 
potential {SV > 0, triangles). The dashed straight line repre- 
sents a dependence expected within the MFA. 



the magnetic impurities. However, as the hole density p 
increases, the data with and without alloy potential tend 
to approach. This is because the dense majority-spin 
holes can embrace more effectively the magnetic impuri- 
ties even without help from the repulsive potential. In 
the other extreme, when the values of p decreases to- 
wards zero, the data with and without alloy potential 
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FIG. 8: Monte Carlo simulations: 
S%.i„{t)/S at different temperatures 



Time evolution of 
1 a Cdo.96Mno.o4Te 



quantum well. The hole concentration is p = 0.4 x 10^^ cm ^ 
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Mn site. 
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approach, too. Here, the Tq magnitude becomes so low 
that the carriers' spins ceases to be entirely saturated in 
the presence of the repulsive potential, which lowers the 
corresponding {\Sz\) value. 



V. INTERPRETATION OF 
MAGNETO-OPTICAL MEASUREMENTS 




T (lO^MCS) 

FIG. 9: Monte Carlo simulations: The magnetization au- 
tocorrelation function G(r) of Mn ions in a Cdo.96Mno.04Te 
quantum well calculated at different temperatures. The con- 
centration of holes p = 0.4 X lO" cm~^ {Nh = 21). Time is 
measured in the Monte Carlo steps per one Mn site. The thin 
solid lines describe exponential decays fitted to the Monte 
Carlo results 



The typical time scale of the observed evolution de- 
pends on many factors including temperature and car- 
rier concentration. To determine this characteristic time 
scale we calculate the time-autocorrelation G(t) of mag- 
netization given 



In this section we present our MC results on Mn mag- 
netization dynamics as a function of temperature and 
carrier concentration. We also discuss whether it is pos- 
sible within the proposed model to estimate the typical 
size of magnetic domains. 



Magnetization dynamics: autocorrelation 
function 



We have traced time dependence of Mn magnetization 
in the absence of an external magnetic field for very long 
runs up to 10^ MC sweeps. We have repeated our sim- 
ulations for different initial configurations including the 
one with all Mn spins aligned along the z-axis. From 
our findings we conclude that except for the first few 
MC sweeps the whole time evolution of magnetization 
does not depend on the initial configuration. The typi- 
cal time-evolution of magnetization at various tempera- 
tures is shown in Fig. [8] For temperatures well below Tq 
the system remains almost unchanged during the whole 
simulation time. For higher temperatures one can ob- 
serve occasional global spin inversions: the system flips 
over between two states of equal energy and spontaneous 
magnetization that differ only in sign. 



g^^-j ^ {m{t)m{t + T)) - {mf 



(12) 



where mit) is the value of magnetization at time t and 
(to) is its average value. The autocorrelation function 
defined by Eq. [12] depends only on the time difference 
T. To calculate this autocorrelation function for up to 
T = 10* MCS computation runs as long as at least 10^ 
MCS are needed. The total computer time depends on 
numerical parameters of which the most important is the 
fc-space cut-off. 

Figure [9]shows the magnetization autocorrelation func- 
tions calculated for the studied system at various tem- 
peratures below and above Tq. Above Tc, i.e., in the 
paramagnetic state, one can observe a very fast decay 
of the autocorrelation function whereas for temperatures 
below Tc correlations persist for times orders of magni- 
tude longer. Within the present MC approach, it is not 
possible to translate numbers of MC steps into physical 
time units in a reliable quantitative way. The qualita- 
tive relations between the time scales of the investigated 
magnetic relaxation processes should nevertheless be re- 
produced in a correct way. Therefore, we undertake only 
a qualitative comparison with the experimental data of 
Secini 
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FIG. 10: Monte Carlo simulations: The magnetization auto- 
correlation function G(t) for selected values of hole concen- 
trations in a Cdo.96Mno.04Te quantum well at T = 0.7 K. 
Time is measured in the Monte Carlo steps per Mn site. 



In order to better understand mechanisms that account 
for magnetization dynamics in a system of competing 
spin-spin AFM and carrier-mediated ferromagnetic inter- 
actions we have also calculated autocorrelation functions 
for a variety of carrier concentration values keeping the 
number of Mn ions unchanged. In Fig. [TO]we display the 
corresponding results calculated for three different carrier 
concentrations at the same temperature T = 0.7 K. As 
one could expect we do observe a significant slowing down 
of magnetization dynamics with increasing the number of 
holes in the QW. Such behavior is understandable within 
the proposed model. Taking into account the magnetic 
anisotropy - all hole-spins align along z-axis - and the 
fact that only the z-component of each Mn-spin couples 
via the p-d exchange to the system of holes, one can see 
that if more holes are present in the system, the value 
of local magnetization and, hence, the barrier height 
for magnetization reversal become higher. Note that 
for the carrier concentrations p = 4.10 



10^° cm-2 and 



p = 8.78 • 10 cm the system is in the ferromagnetic 
state at T = 0.7 K. In contrast, for p = 1.75 • 10^° cm-^ 
the Curie temperature lies below T — 0.7 K (cf. next 
sections), so that the corresponding points describe the 
dynamics of magnetization in a paramagnetic QW. 

A next important aspect one must understand in order 
to interpret properly the experimental findings is the way 
in which the AFM interactions influence the magnetiza- 
tion dynamics. Therefore, we have calculated the magne- 
tization autocorrelation function for the system where all 
short-range AFM interactions are disregarded. In Fig.[TT] 
we compare the obtained magnetic autocorrelation func- 
tions for systems with and without AFM interactions 
calculated at T = 0.7Tc (for each system its own Tc 
as previously determined by MC simulations^ had been 
used). The results presented in Fig. 11 clearly indicate 
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FIG. 11: Monte Carlo simulations: The magnetization au- 
tocorrelation function G(r) of Mn ions in a Cdo.96Mno.04Te 
quantum well at T = Q.lTc- The hole density is p = 
0.4 X lO^'^ cm~^ {Nh — 21). Open circles and full squares: 
antiferromagnetic interactions are neglected {Jij — 0) and 
taken into account (Jij 7^ 0), respectively. Time is measured 
in the Monte Carlo steps per Mn site. 



that short-range AFM interactions may strongly accel- 
erate the decay of the ferromagnetic order. In accord 
with the above observation one can expect that the ef- 
fect of AFM interactions is considerably reduced when 
the layer of Mn ions is thinner than the extend of the 
hole wave function, in other words, when Mn ions are 
concentrated close to the maximum (probability) den- 
sity of holes. To check this we repeated our calculations 
with all Mn ions occupying only the four central layers 
of the investigated QW (note that it is eight layers' thick 
and that for holes the ground state of the inflnite well 
is taken as the "perpendicular" part of their wave func- 
tion, (p{z)). The results of our calculations are shown in 
Fig. [12] One can clearly see that indeed if Mn ions reside 
only in the central layer, the magnetization dynamics re- 
sembles the one for the case when all AFM interactions 
are switched off. We may as well formulate our conclu- 
sion the other way round: more Mn spins close to the 
QW edges, faster relaxation of magnetization. This flnd- 
ing provides a support for theoretical considerations,'^ 
suggesting that magnetization relaxation of bound mag- 
netic polarons in p-type CdTe/(Cd,Mn)Tc QWs occurs 
owing to the AFM coupling to the Mn spins located out- 
side the relevant Bohr radius. 



B. Measuring the domain size: the two-point 
spin-spin connected correlation function 

Let us now look at the spatial distribution of magne- 
tization and, in particular, whether the presence of com- 



11 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



iroooooooooooooooooooooooo 



Mn 



o Mn region reduced 
■ normal QW 



T = 0.7 



Pi 








T (lO^MCS) 



FIG. 12; Monte Carlo simulations: The magnetization au- 
tocorrelation function G(r) of Mn ions in a Cdo.96Mno.04Te 
quantum well at T = 0.7Tc. The hole density is p = 
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R (nm) 



FIG. 13: Monte Carlo simulations: The two-point connected 
correlation function Ci^' (R) of valence band holes at various 
temperatures in a Cdo.geMuo.odTe quantum well. The hole 



0.4 X 10" cm-2 {Nh = 21).0pen circles and full squares: concentration is p = 0.4 x 10 cm {Nh = 21). 
Mn located in the central layers of the QW only and evenly 
distributed within the whole QW, respectively. Time is mea- 
sured in the Monte Carlo steps per Mn site. 



peting short-range antiferromagnetic and long-range fer- 
romagnetic interactions results in the formation of mag- 
netic domains. One of the standard ways to estimate the 
typical domain size is to calculate the spin-spin two-point ^ 
connected correlation function G^{R) given by^ ^ 



-i2/iD\ (m-(r)m(r + R)) — (m)^ 



^5 



GiiR) 



{mn?) — {mY 



(13) 



We have calculated spin-spin correlation functions for 
both: the localized Mn ions and the holes for a number 
of tem pera ture values and various carrier concentrations. 
Figure 



13 



shows the correlation functions (i?) obtained 
for the hole subsystem. We observe long-range correla- 
tions (longer than the size of the simulation box) in the 
ferromagnetic state. For temperatures below Tq the hole 
liquid becomes completely spin-polarized. This has al- 
ready been predicted in the frame of the mean-fi eld ap- 
proximation for the 2D hole gas in a QW.^ In Fig. 14 and 
15 the spin-spin correlation functions for the Mn subsys- 
tem are shown. In Fig. [14] the scale has been stretched 
to make visible how the short-range AFM interactions 
influence the spin-spin correlation function at small dis- 
tances. In this way we can clearly see that for the nearest 
neighbor Mn spins the AFM coupling prevails (the cor- 
relation function becomes negative). At long distances, 
however, below Tc, the long-range ferromagnetic inter- 




FIG. 14: Monte Carlo simulations: The two-point connected 
correlation function g'"c\K) of Mn ions in a Cdo.96Mno.04Te 
quantum well at T = 0.2rc. The hole concentration is p = 
{Nh = 21). 



actions dominate as may be seen in Fig. 15 that shows 
the spin-spin correlation functions on a larger scale. As 
the distance R between the Mn spins becomes large we 
observe (small) positive long-range correlations for tem- 



peratures below Tq. At the same time the correlations 
decrease to zero when temperature increases towards Tq . 
The positive correlations observed in the ferromagnetic 
phase and the lack of correlations in the paramagnetic 
phase may indicate that magnetic domains do exist, but 
their typical size is greater than the simulation box. To 
verify this conjecture further massive simulations would 
be needed. 
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FIG. 15: Monte Carlo Simulations: The long-distance parts of 
two-point connected correlation functions gI?^ (R) of Mn ions 
at various temperatures in a Cdo.96Mno.04Te quantum well. 
The hole concentration is p = 0.4 x lO" cm"^ (Nh = 21). 



VI. CONCLUSIONS 

We have examined magnetization dynamics in a p- 
Cdo.96Mno.04Te quantum well by magneto-optical stud- 
ies and by Monte Carlo simulations taking into ac- 
count the presence of competing short-range antiferro- 
magnetic superexchange and long-range carrier-mediated 
ferromagnetic interactions. Single-particle hole energies 
and eigenfunctions for particular Mn spin configurations 
have been updated numerically after each Monte Carlo 
sweep. In addition to the p-d exchange interaction, spin- 
independent alloy disorder has been taken into account 
and found to be important. In particular, a sufficiently 
strong attractive potential introduced by the magnetic 
ion leads to hole localization and to the corresponding 
disappearance of a long-range order mediated by the car- 
riers. Also a strong repulsive potential, by reducing the 
magnitude of the carrier wave function at the magnetic 
ion, may diminish the magnitude of Tq. This effect is 
particularly strong in the case of a quantum well, where 
the direction of the hole spin is fixed rather by the spin- 
orbit interaction than by the p-d coupling, so that the 



repulsive potential may not be compensated by the p-d 
interaction. 

Both experimental and theoretical results demonstrate 
that the transition from the paramagnetic to ferromag- 
netic phase increases the magnetization relaxation time. 
In particular, the relaxation time determined from trac- 
ing the time evolution of magnetization after a pulse of 
the magnetic field is shorter than 20 ns in the paramag- 
netic state. Lowering temperature below Tq results in 
an increase of the relaxation time up to only 2 fis. Our 
numerical results indicate that antiferromagnetic inter- 
actions between Mn spins account for this relatively fast 
magnetization fluctuations and, thus, might be responsi- 
ble for the absence of spontaneous magnetization below 
Tc, as found in static measurements. While our Monte 
Carlo simulations of the two-point spin-spin connected 
correlation function show a huge influence of the short- 
range antiferromagnetic interactions on the magnetiza- 
tion dynamics, they do not reveal the formation of do- 
mains that would be smaller than our simulation cell, 
L = 226 nm. 

At the same time, our Monte Carlo results reveal that 
magnetization relaxation in the ferromagnetic phase pro- 
ceeds primarily due to antiferromagnetic couplings to Mn 
spins residing at the quantum well boundary, as they are 
weakly polarized by the carriers. Accordingly, magneti- 
zation dynamics would be much slow down if the width 
of the region containing Mn spins were been narrower 
than the extent of the carrier wave function. We predict, 
therefore, that spontaneous magnetization and the asso- 
ciated magnetic hysteresis could be observed in quantum 
wells containing the Mn layer only in its center. 
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